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Abstract 

Objective: The study of obesity-related metabolic syndrome or Type 2 diabetes (T2D) in children is particularly difficult 
because of fear of needles. We tested a non-invasive approach to study inflammatory parameters in an at-risl< population of 
children to provide proof-of-principle for future investigations of vulnerable subjects. 

Design and Methods: We evaluated metabolic differences in 744, 1 1-year old children selected from underweight, normal 
healthy weight, overweight and obese categories by analyzing fasting saliva samples for 20 biomarkers. Saliva supernatants 
were obtained following centrifugation and used for analyses. 

/feiu/ts; Salivary C-reactive protein (CRP) was 6 times higher, salivary insulin and leptin were 3 times higher, and adiponectin 
was 30% lower in obese children compared to healthy normal weight children (all P<0.0001). Categorical analysis 
suggested that there might be three types of obesity in children. Distinctly inflammatory characteristics appeared in 76% of 
obese children while in 13%, salivary insulin was high but not associated with inflammatory mediators. The remaining 1 1% 
of obese children had high insulin and reduced adiponectin. Forty percent of the non-obese children were found in groups 
which, based on biomarker characteristics, may be at risk for becoming obese. 

£ro/7c/t/5/b/75; Significantly altered levels of salivary biomarkers in obese children from a high-risk population, suggest the 
potential for developing non-invasive screening procedures to identify T2D-vulnerable individuals and a means to test 
preventative strategies. 
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introduction 

The prevalence of pediatric obesity lias increased worldwide in 
recent years and raised urgent concern about expected metabolic 
dysregulation and serious comorbidities, such as T2D, likely to 
arise among obese children as they age into adulthood [1]. The 
relationship is unknown among adult-onset, obesity-associated 
pathologies and exposure of these individuals to diet-induced 
obesity while still children. Furthermore, it is not known if at-risk 
cohorts of pediatric subjects can be identified by non-invasive 
measures early enough during development of pediatric obesity for 
effective interventions to reduce risk of T2D. The importance of 
the answers to these questions is widely recognized, thus countries 
with significant prevalence of pediatric obesity are devoting 
resources to understand the scope of this problem and devise 
appropriate public health recommendations. 



In this research, we describe an approach that uses saliva 
samples as a basis for diagnosis. The selection of biomarkers for 
evaluation in this study of obese children was focused on 
adipokines (adipocytokines), because concentrations of these 
factors partly define metabolic health. A subset of these proteins 
is of particular interest for immunometabolism, including 
adiponectin, leptin and resistin, which are adipocyte-derived 
hormones secreted in saliva [2] , salivary [3] monocyte chemotactic 
protein- 1 (MCP-1) and salivary [4,5] tumor necrosis factor-alpha 
(TNF-ot). Interleukins (IL) comprise a second major class of 
investigated biomarkers, each of which has also been measured in 
saliva, but not in an obese adolescent population. Specifically, IL- 
1 P is associated with periodontal inflammation [6] ; IL-6 has been 
measured in periodontitis patients; IL-4, IL- 1 0, IL- 1 2 and IL- 1 7 
are related to Sjogren's syndrome [7]; IL-10 is reduced in 
periodontitis patients [5]; IL-8 is related to dental caries in 
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adolescents [4] and oral cancer [8]; IL-13 is identified in the 
sputum of asthmatics [9]; IL-17 is lower in patients with 
periodontal disease [10]; and interferon y (IFN-y) is higher in 
the saliva of control subjects without Sjogren's syndrome [7]. 
Thus, adipokines and cytokines from sources other than peripheral 
blood have been assayed and associated with disease states. Both 
salivary immunoreactive insulin [1 1] and salivary ghrelin [1 2] 
have been associated with T2D and obesity, and were included in 
our analysis. Two oral inflammatory disease markers were also 
included, matrix metalloproteinase-9 (MMP-9) and myeloperox- 
idase (MPO). MMP-9 is a protease often found elevated in patients 
with periodontal disease [13] and oral cancer [14]; whereas MPO, 
a peroxidase often used as a measure of neutrophil degranulation, 
is elevated in diabetic patients [15]. Vascular endothelial growth 
factor (VEGF) was included because it appears in saliva of diabetic 
women [16]. C-reactive protein (CRP) was included because it is 
associated with inflammation and cardiovascular disease [17]. 
Collectively, the panel defines inflammatory and metabolic 
processes important for metabolic health in obesity. 

Blood elements that partition into saliva have attracted interest 
[18]. Saliva has major advantages over peripheral blood for the 
study of disease development in childhood because sample 
collection is non-invasive. Here, we collected 8,319 saliva samples 
from 10-12-year old Kuwaiti children, a high risk population for 
obesity and T2D development [19]. We measured salivary 
concentrations of 20 hormones and cytokines in a random subset 
of 744 subjects. The well-established relationship between insulin- 
resistant obesity and systemic elevation of pro-inflammatory 
cytokines [20] supports our hypothesis that the saliva of obese 
children shows informative imbalances in critical immunometa- 
bolic factors [21]. The results of this study provide insight into the 
early development of metabolic disease in children and establish 
that non-invasive methods are robust and useful for data collection 
in studies of vulnerable subjects. 

Materials and Methods 

Subject Selection 

The validation study for US subjects was approved by the 
Forsyth Institutional Review Board in Cambridge MA, U.S.A. 
The study of Kuwaiti subjects (4''^ or 5* grades) was approved by 
the Dasman Diabetes Institute Ethical Review Committee in 
Kuwait. Arabic language informed consent was signed by parents/ 
guardians in advance. Subject assent was obtained the day of the 
visit. 

Data Collection 

Data and saliva samples were collected from 8,319 subjects 
during 182 visits to 39 Kuwaiti schools between October 2, 201 1 
and May 15, 2012. Subject identification, height, weight, blood 
pressure, food preferences, oral examination, fitness and sleep 
parameters were entered into a programmed iPad (Apple, 
Cupertino, CA) system for internet transfer. Fitness was measured 
by heart rate elevation (beats/minute) following a standard 3- 
minute exercise [22]. Body weight categories were defined using a 
Body Mass Index (BMI) z-score [23]. Blood and saliva samples 
were collected for insulin analysis from February 23-September 
23, 201 1 from 53 US subjects living in Maine and Massachusetts. 
Otherwise, the US and Kuwait protocols were identical. 

Saliva Collection 

After a 15 ml water rinse of the mouth, whole saliva was 
collected by expectoration of approximately 3 ml into a sterile 
tube kept on ice. Samples were collected starting at 8:30 AJvI 
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Table 2. Concentration of insulin, C-reactive protein (CRP), adiponectin and leptin in the saliva supernatant of 744 children by 
body weight category and sex. 





Probe 


insulin 


CRP 


Adiponectin 


Leptin 




Parameter 


(median, Intraquartile 
range) 


(median, intraquartile 
range) 


(median, intraquartile range) 


(median, intraquartile 
range) 




Units 


pg/ml saliva 


pg/ml saliva 


pg/ml saliva 


pg/ml saliva 


Underweight 
(1 86) 


M(64) 


39.17, 48.63 


56.64, 92.32 


4,421, 6,424 


0.01, 5.31 




F(122) 


34.30, 38.02 


61.99, 182.74 


5,060, 4,573 


0.01, 3.22 


Normal (186) 


IVI(93) 


39.39, 45.38 


73.01, 153.75 


4,220, 5,303 


1.06, 4.77 




F(93) 


44.70, 54.38 


77.15, 186.95 


3,994, 5,052 


0.63, 4.61 


Overweight 
(186) 


M(93) 


80.39, 88.74 


177.46, 311.93 


2,402, 3,785 


1 .06, 3.26 




F(93) 


76.25, 87.13 


281.39, 516.54 


3,322, 3,693 


2.41, 5.03 


Obese (186) 


IVI(93) 


112.98, 125.09 


429.44, 668.52 


2,548, 2,779 


3.16, 6.40 




F(93) 


143.50, 150.24 


443.13, 1033.29 


3,062, 3,752 


3.70, 6.41 


Wllcoxon 
regression p 


Age 


0.107 


0.083 


0.373 


1.000 




Sex 


0.337 


0.028 


0.471 


1.000 




Overweight 


< 0.0001 


<0.0001 


0.001 


0.408 




Obese 


< 0.0001 


< 0.0001 


< 0.0001 


<0.0001 




Underweight 


0.157 


0.266 


0.142 


1.000 


Summary statistics are median, interquartile range (N subjects). Values recorded for sixteen additional probes tested are in Table S3 in File SI. Probability levels for age, 
sex, overweight, obese and underweight were computed by Wilcoxon regression relative to normal healthy weight children. 



doi:l 0.1 371 /journal.pone.0098799.t002 



before the subjects had breakfast. The time of start and end of 
collection was recorded for all children. Iced samples were 
immediately transported to the Dasman Institute where they were 



weighed to determine volume collected and centrifuged to remove 
particulate debris and exfoliated mucosal cells at 2,800 RPM for 
20 minutes at 4°C. Supernatant aliquots were transferred to 



Insulin 
pg/ml 



200 
150 
100 
50 
0 



Leptin 
pg/ml 




1200r 



CRP 

pg/ml 



800 



400 



Under- Normal Over- Obese 
weight weight weight 




7000 



Adiponectin 
pg/ml 



5000 



Under- Normal Over- Obese 
weight weight weight 



3000 



lOOOL 




Under- Normal Over- Obese 
weight weight weight 




Under- Normal Over- Obese 
weight weight weight 



Significantly different from normal weight (p < 0.001) 



Figure 1. The concentration of blomarkers in saliva supernatant that were significantly different in obese children compared to 
normal weight children. Insulin, CRP and leptin significantly increased with obesity whereas adiponectin decreased. Central values represent the 
median concentration and whiskers represent the interquartile range (4-75* percentile, -25* percentile). 
doi:1 0.1 371/journal.pone.0098799.g001 
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screw-cap 2D barcoded storage tubes (Thermo Scientific) read by 
a barcode reader (Thermo Scientific VisionMate ST). Tlie 
barcode was captured with subject number to spreadsheet. The 
sample vials were sealed by a torque-controUed tube capper 
(Thermo Scientific 8-Channel Screw Cap Tube Capper), placed in 
a 96-vial rack (Thermo Scientific Latch Rack) and fi-ozen at — 
80°C. Racks were air-transferred fi'om Kuwait under temperature- 
monitored dry ice (Biocair, Boston MA) to the Forsyth Institute 
and maintained at — 80°C until assay (average time to 
assay = 0.88±0.06 y). 

From the 8,319 samples, 744 samples were selected for assay 
such that 93 were from subjects of each of the defined body weight 
categories described in the data analysis section for each sex. 
Random selection was achieved by assigning a random number 
between 1 and x = integer (total number in each group/93) to each 
subject, sorting for all assigned the number 1 and selecting the first 
93. All body weight categories included 186 subjects evenly 
distributed by sex except for underweight subjects, where, due to a 
shortage of male underweight subjects, 64 male and 122 female 
subjects were selected. 

Multiplex Analysis of Salivary Markers 

AH assays (744 saliva samples for 20 biomarkers) were 
performed on 25 \ll of saliva supernatant using four multiplex 
magnetic bead panels on a Luminex 200 platform (Luminex, 
Austin, TX). Results were evaluated using Bio-Plex Manager 
(Version 5.0; Bio-Rad, Hercules, CA). IFN-y, IL-10, interleukin- 
12p70 (IL-12P70), IL-13, IL-17A, IL-ip, IL-4, IL-6, IL-8, MCP-1, 
TNF-a and VEGF were measured by a 12-plex human cytokine/ 
chemokine panel with no dilution (MiUipore cat ^HCYTOMAG- 
60K; lot T^2055690). Metabolic hormones ghrelin, insulin and 
leptin were measured by a 3-plex human metabolic hormone 
panel with no sample dilution (MUlipore cat #HMHMAG-34AK; 
lot #2055724). MPO and MMP-9 were measured by a 2-plex 
human cardiovascular disease panel with 1:2 sample dilution 



200 



150 



Saliva Protein 
Concentration 
(mg/DL) 



100 



Protein = -0.57 x Salivary flow rate + 66.7 

• * r=0.30 

p=0.00001 




20 40 60 80 
Salivary flow rate 
(ml/hour) 



100 



Figure 2. Total salivary protein concentration and salivary flow 
rate. In the subset of 213 samples from children measured for total 
salivary protein, protein content significantly declined with increased 
salivary flow rate. 

doi:10.1371/journal.pone.0098799.g002 

(Millipore cat #HCVD1-67AK; lot #2055723). Adiponectin 
(total), CRP and resistin were measured in a 3-plex human obesity 
panel with 1:2 sample dilution (R&D cat #LOB000, LOB1065, 
LOB 1707, LOB 1359; lot #300710). The panel used to evaluate 
insulin in blood and saliva of US subjects was a human metabolic 
hormone panel single plex for insulin (Millipore cat #HMHMAG- 
34K-01EMD) with no dilution. All assays were performed 
following manufacturers' protocol, with the exception of an 
additional 3 standards to increase the range of detection [24]. 



Table 3. Concentration of insulin, CRP, adiponectin and leptin normalized by total protein of a subset of 213 randomly selected 
children. 





Probe 


Insulin 


CRP 


Adiponectin 


Leptin 






(median, intraquartile 


(median, intraquartile 


(median. 


(median, intraquartile 




Parameter 


range) 


range) 


intraquartile range) 


range) 




Units 


pg/ml saliva 


pg/ml saliva 


pg/ml saliva 


pg/ml saliva 


Underweight (27) 


M(12) 


39.11, 84.01 


61.66, 143.95 


13, 363, 13, 812 


0.02, 4.47 




F(15) 


87.84, 83.83 


117.67, 217.16 


8, 168, 6, 591 


0.03, 4.66 


Normal (71) 


M(36) 


86.45, 1 33.87 


142.82, 236.09 


8, 919, 7, 590 


1.04, 7.65 




F(35) 


53.56, 83.21 


180.74, 327.00 


11, 426, 10, 123 


0.67, 8.59 


Overweight (57) 


M(19) 


192.29, 208.82 


430.76, 625.88 


7, 152, 5, 916 


5.14, 5.82 




F(38) 


168.09, 257.48 


607.19, 787.82 


9, 184, 6, 897 


3.38, 8.71 


Obese (58) 


M(27) 


211.83, 203.68 


1, 038.31, 1, 484.94 


6, 337, 3, 550 


3.39, 12.44 




F(31) 


188.95, 236.98 


828.61, 1, 560.40 


7, 740, 5, 233 


4.77, 11.75 




Age 


0.047 (0.038) 


0.825 (0.947) 


0.080 (0.097) 


0.970 (1.000) 


Wilcoxon regression p 


Sex 


0.452 (0.619) 


0.564 (0.253) 


0.181 (0.329) 


0.970 (1.000) 


mg protein (ml saliva} 


Overweight 


0.262 (0.448) 


0.132 (0.266) 


0.560 (0.911) 


0.018 (0.020) 




Obese 


<0.0001 (<0.0001) 


<0.0001 (<0.0001) 


0.082 (0.030) 


0.01 7 (0.020) 




Underweight 


<0.0001 (<0.0001) 


<0.0001 (<0.0001) 


0.028 (0.017) 


0.003 (0.012) 



Probability levels for protein-normalized saliva samples for age, sex, overweight, obese and underweight were computed by Wilcoxon regression relative to normal 
healthy weight children are listed with probability levels of the subset based on salivary concentration in parentheses. 
doi:l 0.1 371 /journal.pone.0098799.t003 
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43.60 




Median 


22.13 
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Figure 3. Salivary protein concentration and salivary flow rate of children by body weight category. Neither protein content nor salivary 
flow rate differed significantly in this subset of 213 samples between children in each body weight category. 
doi:1 0.1 371 /journal.pone.0098799.g003 



Saliva supernatants were thawed at 4°C overnight and kept on 
ice throughout the assay procedures. Briefly, all kit components 
were brought to room temperature and reagents were prepared as 
instructed (wash bufiers, beads, standards, etc.). Assay plates (96- 
well) with assay buffer, standards, samples and beads were covered 
and incubated on plate shaker (500 rpm) for either 3 hours at 
room temperature for the obesity panel, or overnight at 4°C. 
Plates were washed and detection antibody cocktail added. Plates 
were covered and incubated at room temperature for 1 hour on a 
plate shaker. Streptavidin-phycoerythrin fluorescent reporter was 
added to all wells, and the covered plate incubated for 30 minutes 
at room temperature on a shaker. Plates were washed three times 



and beads resuspended in wash buffer, placed on shaker for 
5 minutes, read and evaluated. 

Salivary Protein Assay 

A subset of 2 1 3 saliva supernatants (approximately 54 from each 
body weight category) were randomly selected and assayed for 
total protein. The assay was by a commercial kit (Pierce 660 nm 
protein assay, Pierce Biotechnology Rockford, IL) which has a 
manufacturer's stated working range of 50-2,000 ]Xg/ ml. Values 
were determined on 1 0 |J,1 saliva supernatant samples using bovine 
serum albumin as a standard implemented on a robotic chemical 




(pg/ml) 



Figure 4. A: Categorical decision tree describing identification of children that are obese or not obese. 76% of the obese children 
were identified by the predictor variable salivary CRP >219 pg/ml. Of the lower salivary CRP saliva samples, 13% of the obese children were 
identified by the predictor variable insulin >128 pg/ml and 11% had insulin £128 pg/ml. B: Insulin concentration in fasting saliva supernatant and 
plasma of 53 adolescent donors. These data were fitted after log transformation to obtain the equation Ln(Plasma) = 0.85*Ln(Saliva) +1.84 with 
r^ = 0.67. Based on this approximation, the predictor variable of 128 pg/ml saliva insulin would be approximately 67 pmoles/L of plasma insulin (or 
11 |J.U/ml using the conversion factor 1|j.U/ml = 6.00 pmol/L). 
doi:1 0.1 371/journal.pone.0098799.g004 
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Table 4. Median values for measures and biomarkers of groups defined by CRP and insulin concentration in saliva supernatant. 







Non~inf iBrnmcd 








Inf iBrnmsd 






CRP S219 








CRP >219 






Insulin £128 




Insulin >128 








Clinical Measure 


OH 


NH 


Ol 


Nl 


OC 


NC 


BMI (Kg/M^) 


26A*** 


16.2 


29.8*** 


21 .4*** 


28.2*** 


21.1*** 


Waist circumference (cm) 


78.7*** 


58.4 


86.4*** 


68.6*** 


83.8*** 


66.0*** 


Systolic BP (mmHg) 


121*** 


101 


122*** 


113** 


1 20*** 


1 09*** 


Diastolic BP (mmHg) 


78* 


70 


84*** 


75* 


79*** 


72.0 


Fitness (beats/min) 


34.8** 


19.5 


29.5 


20.0 


35.5*** 


23.0 


N 


20 


337 


24 


44 


142 


177 


% of total 


2.7 


45.3 


3.2 


5.9 


19.1 


23.8 


% of obese 


10.8 




12.9 




76.3 






CRP <219 








CRP >219 






Insulin £128 




Insulin >128 








Biomarker 


OH 


NH 


Ol 


Nl 


OC 


NC 


CRP (pg/ml) 


62.3 


54.1 


126.6** 


67.1 


604.6*** 


521.0*** 


Insulin (pg/ml) 


74.5*** 


41.8 


186.8*** 


1 77.4*** 


121.5*** 


54.3*** 


IL-6 (pg/ml) 


3.8 


6.2 


3.4* 


4.6 


8.4*** 


10.7*** 


Leptin (pg/ml) 


1.3 


0.6 


5.7** 


0.5 


3.3*** 


1.4* 


IL-10 (pg/ml) 


2.8*** 


4.3 


2.7* 


2.8* 


5.5* 


6.1*** 


Adiponectin (pg/ml) 


1201*** 


3582 


1872** 


1788** 


3134 


4816** 


Resistin (pg/ml) 


1073* 


1555 


1116 


948* 


1812 


2296*** 


IL-ip (pg/ml) 


24.5 


41.9 


16.4* 


24.7* 


47.7 


77.6*** 


MMP-9 (pg/ml) 


8.9* 


17.1 


8.3* 


12.8 


22.7* 


27.5*** 


IL-8 (pg/ml) 


260.5* 


440.8 


243.4* 


259.4* 


476.7 


644.8** 


TNF-a (pg/ml) 


4.3* 


5.8 


1.8** 


3.1** 


6.6 


7.9** 


VEGF (pg/ml) 


206.5* 


265.3 


308.4 


348.8** 


290.5* 


317.7* 


IL-12P70 (pg/ml) 


13.2 


14.7 


8.2 


9.3 


9.9 


17.5* 


IL-1 7A (pg/ml) 


1.6 


1.0 


0* 


0.2* 


0.7 


1.5* 


MCP-1 (pg/ml) 


174.1* 


236.8 


282.3 


308.3* 


267.0* 


259.8 


Ghrelin (pg/ml) 


1.9 


1.1 


1.9 


1.1 


1.1 


1.4 


IFNy (pg/ml) 


3.0 


2.2 


1.0* 


1.3 


1.9 


2.6 


IL-1 3 (pg/ml) 


6.1 


6.9 


3.0* 


4.5* 


6.7 


7.1 


MPO (pg/ml) 


263 


30.7 


25.5 


32.5 


30.5 


31.7 


IL-4 (pg/ml) 


7.7 


7.5 


5.2 


5.8 


7.2 


8.8 



Categorical comparisons were by the Mann-Whitney U test. A Bonferroni correction of 0.0004 (for overall p<0.05) was applied to control the familywise error rate 

associated with 125 test procedures (in boldface). 

***p<0.0004, 

**p<0.001, 

*p<0.05. 

doi:1 0.1 371 /journal.pone.0098799.t004 



analysis unit (EVO 100 with a Sunrise spectrophotometer set at 
660 nm, Tecan US, Morrisvile NC). 

Data Analysis 

Body weight categories were defined on BMI z-score based on 
pubhshed growth charts [23]. By this criterion, obese was &95* 
percentile, overweight was &85*— 95* percentile, normal healthy 
weight was S:5th and 85* percentile and underweight was <5* 
percentile [25]. Salivary flow rate was computed from the start 



time, stop time and tarred salivary weight. To account for the non- 
normality of multiprobe data, a multivariate rank-based WUcoxon 
regression method [26] was applied. The relationship between the 
probe concentration and body weight categories were adjusted for 
age and sex. Values for age, BMI, waist circumference (WC), 
systolic and diastolic blood pressure, and fitness were analyzed for 
significance between body weight categories by a Kruskal-WaUis 
rank sum test followed by pairwise comparisons using the 
Wilcoxon rank sum test. Predictor variables, diagnostic sensitivity 
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and specificity were determined by CART software (Salford 
Systems, San Diego, CA). Log transformation of both insulin 
measurements preceded their analysis. 

Results 

Saliva samples (744) were randomly selected from the 8,319 
total to provide 186 for each category of underweight, normal 
healthy weight, overweight and obese. Characteristics of these 
groups are evaluated in Table 1 . No significant differences in age 
between groups were noted. BMI, WC, blood pressure and fitness 
exhibited significant incremental increases as the body weight 
category moved to obese. In every body weight category, BMI and 
WC significantiy differed from each other. Systolic and diastolic 
blood pressures were significantly different only between normal 
and overweight or obese subjects. Obese subjects had significandy 
lower fitness as measured by exercise stimulated heart rate 
eli'\'atiori than any other category. Obese subjects had 64% 
higher BMI, 40% greater WC, 17% higher systolic and diastolic 
blood pressure and 50% greater exercise elevation of heart rate as 
a measure of unfitness compared to normal weight subjects. Table 
SI in File SI shows similar statistics from 53 US subjects used to 
measure the relative concentration of insulin in saliva and plasma. 

The assay performance is summarized in Table S2 in File SI. 
The median concentration of 1 7 biomarkers was greater than the 
manufacturer's stated assay sensitivity. Median values fell below 
the assay sensitivity for IL-10, leptin and ghrelin. The analysis 
software allowed extrapolation beyond the lowest standard, 
provided the value did not fail on or below the blank. Using this 
feature, although the concentration cannot be said to have been 
accurately measured, the assay provided a means to order low 
concentrations for non-parametric rank analysis. All probes 
provided non-zero values for 95 to 100% of the samples tested 
except for leptin; 37.4% of the leptin determinations were too low 
to measure. In the leptin analysis, however, the percentage of 
measurable samples increased with increasing obesity so that 
analysis by non-parametric rank could be performed without 
introducing bias. 

The assay results are summarized in Figure 1 and Table 2. Of 
the 20 probes tested, significant differences between obese and 
normal weight subjects occurred in the concentration of insulin, 
CRP, adiponectin and leptin. CRP and insulin were significantiy 
elevated in obese relative to normal subjects. Salivary insulin levels 
in obese subjects (median =127 pg/ml) were almost 3 times that of 
normal weight subjects (median = 44 pg/ ml) and salivary CRP of 
obese subjects (median = 435 pg/ml) was almost 6 times that of 
normal weight subjects (median =76 pg/ml). The estimated 
concentration of salivary leptin in obese subjects (median = 3.3 
pg/ml) was 3 times that of normal weight subjects (1 pg/ml). 
Salivary adiponectin decreased by approximately 30% with 
increasing obesity from that of normal healthy weight subjects 
(median 4,083 pg/ml) to overweight and obese subjects (medi- 
an =2981 and 2798 pg/ml respectively). Differences between 
normal and obese subjects for insulin, CRP, adiponectin and 
leptin were all highly significant (p<0.0001). There were no 
significant differences between biomarker concentrations of 
underweight and normal healthy weight subjects. Salivary levels 
of the remaining 16 probes (Figure SI, Table S3 in File SI) did not 
significantiy change when analyzed by body weight categories. 

Evaluation of saliva supematants normalized by salivary total 
protein content (Table 3), indicated that significant levels of 
biomarkers obtained were similar to those measured using salivary 
concentration. Salivary protein concentration significantiy de- 
creased with increasing salivary flow rate (Figure 2). Neither 



safivary protein concentration nor salivary flow rate varied 
significantiy between body weight categories (Figure 3). 

The diagnostic implications of these findings were evaluated by 
classification tree topology (Figure 4A). Of the obese subjects, 
76.3% had >219 pg/ml CRP in their saliva, indicating that the 
inflammatory state was the most common form of obesity. Of the 
obese subjects with lower levels of CRP (<219 pg/ml), 13% had 
high insulin (>128 pg/ml) and 11% had low insulin. The obese 
group with CRP <219 and insulin <128 would be predicted to be 
'metabolically healthy'. Using only salivary CRP and insulin as 
predictors, the overall diagnostic sensitivity for identifying obesity 
by biomarker status was 89% and specificity was 61%. 

We used Luminex analysis to evaluate correspondence between 
saliva and plasma concentration of insulin in a group of 53 
American subjects similar in age to the Kuwaiti subjects 
(Figure 4B). We found that the concentration of immunoreactive 
saliva insulin is approximately half of the plasma concentration 
and estimates of plasma insulin with high correlation (r^ = 0.67). 
By this approximation, the predictor variable of 128 pg/ml 
salivary insulin is equivalent to 68 pmoles/L = 1 1 |i,U/ ml in 
plasma, a value higher than that reported for normal fasting 
subjects 38-46 pmoles/L = 6.3-7.6 |lU/ml [27]. 

Investigation of the properties of each of the six groups 
described by this categorization tree revealed significant differ- 
ences relative to the normal healthy (NH) group in both clinical 
measures and biomarker concentrations (Talilc 4). Systolic blood 
pressure, BMI, WC and insulin were significantiy elevated in all 
categories except for the non-obese healthy group (NH). Diastofic 
blood pressure was significantiy elevated in the obese with high 
insulin (OI) and obese with high CRP (OC) groups. Low fitness 
was significantly increased in the OC group. CRP and IL-6 were 
significantiy elevated in both high CRP groups (OC and NC) and 
leptin was significantiy elevated in the obese high CRP group 
(OC). The non-obese, high CRP group (NC) exhibited elevation of 
multiple biomarkers (IL-10, resistin, IL-lfi and MMP-9). The 
obese but healthy group (OH) exhibited significantiy lower 
concentrations of IL-10 and adiponectin. 

Discussion 

Many blood factors also occur in saliva and their levels often 
correlate. In our study, we surveyed 20 possible biomarkers related 
to obesity and found four that exhibit significant change with 
increasing body weight in a pediatric population. These data 
suggest that saliva could be a useful blood surrogate for the study 
of metabolic complications of obesity in children, where repeated 
blood sampling can be both traumatic and difficult. 

Elevated plasma insulin is characteristic of T2D and also is 
proportional to body fat content [28] . In agreement with our data, 
investigators have noted that salivary insulin levels in obese 
children may be four or five times higher than in normal healthy 
weight children [29] and that salivary insulin correlates well with 
plasma levels at only slightly lower concentration [30]. Salivary 
insulin is linearly related to plasma insulin during the glucose 
tolerance test [1 1] and correlates well with plasma concentration 
after insulin injection [31]. Plasma insulin decreases in parallel 
with weight loss in obese children [32]. 

Salivary CRP concentration has also been found to correlate 
well with serum concentrations [24] . Association between salivary 
CRP and obesity has been reported in black South African 
children [33]. Human studies associate high levels of CRP with 
metabolic syndrome and T2D [34]. Reduced adiponectin is an 
independent risk factor for metabolic syndrome [35]. Both salivary 
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adiponectin [24] and leptin [2] concentrations are highly 
correlated with plasma concentrations. 

Since it is well-known that the concentration of proteins can 

var^' significantly in response to stimulation, total protein is 
sometimes used to normalize concentration of specific salivary 
analytcs [36]. The total protein content of salivary samples was 
reduced as salivary flow rate increased (Figure 2), however, neither 
total salivary protein nor sahvary flow rate were significantly 
different between body weight categories (Figure 3) indicating that 
obesity did not affect salivary protein levels or salivary flow rate in 
this population. Investigation of a protein-normalized sample 
subset indicated that this procedure did not appreciably modify 
the diagnostic potential of salivary insulin, CRP, leptin or 
adiponectin in the samples collected from adolescent Kuwaiti 
children (Table 3). 

The first level of the classification tree (CRP >219 pg/ml. 
Figure 4A), is associated with multiple inflammatory mediators 
(Table 4). The obese high salivary CRP (OC), represented the 
largest group of obese subjects (76.3%). The major differences 
among groups were elevated diastoUc blood pressure, significandy 
reduced fitness and elevated leptin in the obese (OC) group. The 
non-obese high salivary CRP group (NC) had high levels of CRP, 
insulin, IL-6, IL-10, resistin, IL-ip and MMP-9 compared to the 
normal healthy group (NH). 

The second level of the classification trcx; (CRP <219, insulin > 
128 pg/ml) also identified two groups of subjects (OI and NI), but 
these did not have elevated salivary inflammatory mediators but 
did have elevated insulin le\-fls. Th(^ principal differences were that 
obese subjects had high systolic and diastolic blood pressures. 
Subjects with both low insulin and low CRP included a normal 
healthy group (NH) and an obese, 'healthy' group (OH). In this 
case, obese subjects were associated with elevated insulin, reduced 
adiponectin and reduced IL-10. 

Comparison of these observations with those based on 
immunohistochemistry of obese adult adipose tissue [37] strongly 
suggests that the group designated as CLS+ (macrophage crown- 
like structure positive), found in 72% of obese adults, may coincide 
with the subjects designated as CRP >219 in our studies (76.3% of 
obese subjects). Our studies suggest, however, that non-obese 
cohorts with comparable biomarker levels could be at risk for 
development of immunometabolic complications of obesity. 

The constellation of biomarkers that appear in saliva are likely a 
reflection of underlying pathology [38-40] . High concentrations of 
salivary CRP, myoglobin and MPO, for example appear after 
myocardial infarction [41]. Periodontal disease characteristically 
exhibits high levels of salivary lL-8 [42], MMP-9 and IL-ip [43]. 
Inflammatory bowel disease is associated with elevated salivary IL- 
6 and CRP [44] . We can now conclude that adolescent obesity is 
most commonly associated with high levels of salivary insulin and 
CRP. 

Conclusions 

We have identified four salivary biomarkers in 10-year old 
subjects that change significandy with increasing obesity; insulin, 
CRP, adiponectin and leptin. The results of this study suggest that 
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